Duality with real-space renormalization and its application to bond percolation 
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We obtain the exact solution of the bond-percolation thresholds with inhomogenous probabilities 
on the square lattice. Our method is based on the duality analysis with real-space renormaliza- 
tion, which is a profound technique invented in the spin-glass theory. Our formulation is a more 
straightforward way compared to the very recent study on the same problem [R. M. Ziff , et. ai, J. 
Phys. A: Math. Theor. 45 (2012) 494005]. The resultant generic formulas from our derivation can 
give several estimations for the bond-percolation thresholds on other lattices rather than the square 
lattice. 



I. INTRODUCTION 

Forest fire happens suddenly and spreads out rapidly. 
In order to save the forest itself, living animals, humans 
and their community there, it is important to resolve 
a naive question: how can we prevent the fire spread 
through the whole system? In the present study we take 
an associated mathematical problem, namely percolation 
[l| . The percolation is a very simple but ubiquitous prob- 
lem, which is closely related to the phenomenon involved 
in the formation of long-range connectivity in systems, 
as well as forest fire as exemplified above. For instance, 
it provides rich comprehensions for numerous practical 
issues including conductivity in composite materials, in- 
fectious disease, flow through porous media, and poly- 
merization. In the present study we restrict ourselves 
to the case of the bond percolation problem, where each 
bond to connect both ends on the system is selected in 
a stochastic manner. The bond percolation is a typi- 
cal instance of the cooperative phenomena, with which 
highly skillful techniques are essential to deal. Neverthe- 
less very simple formulas have been expected to hold for 
the bond-percolation thresholds at which giant clusters 
over the whole system appear. The key is a particular 
symmetry embedded in the system, namely the duality. 

In classical spin models, the duality is known to be 
a hidden symmetry between the partition functions in 
low and high temperatures. This symmetry allows us 
to identify the locations of the critical points for various 
spin models such as the Ising and Potts models 0, Q. 
In the present study we employ the duality in order to 
assess the bond-percolation threshold, since g-state Potts 
model can be mapped to the bond-percolation problem 
in the limit of q — > 1 0, Q ■ The special symmetry of the 
square lattice, namely self-duality, yields the exact solu- 
tion of the bond-percolation threshold in the case with a 
homogenous probability on each bond. Even for the case 
without self-duality, we can perform the duality analy- 
sis to obtain the bond-percolation thresholds in several 
cases in conjunction with another technique, namely the 
star-triangle transformation [ij ■ 

In the present study, we generalize the star-triangle 
transformation to the case on the square lattice. We ap- 
ply the generalized technique, namely the duality analy- 



sis with real-space renormalization, to the inhomogenous 
case on the square lattice. The resultant equation to pro- 
vide the bond-percolation thresholds coincides with that 
proposed by Wu @ . Very recent work performed by Ziff, 
et. al. has proved its validity by combination of the 
several profound results Q. Our technique provides a 
more straightforward way to derive the exact formula on 
the critical manifolds of the bond-percolation thresholds 
without any other ingredients to support our analysis. 
Moreover, we give explicit forms of several generic for- 
mulas depending on the structure of the unit cell form- 
ing the lattice. The basis of our technique comes from 
the different stream of study on random spin systems, in 
particular spin glasses. The straightforward rederivation 
of the existing equalities in the different context implies 
existence of close connection between different realms, 
bond-percolation problems and spin glasses. 

The paper is organized as follows. In the next sec- 
tion, we review the conventional duality and the star- 
triangle transformation for convenience. The third sec- 
tion demonstrates the duality with real-space renormal- 
ization to the inhomogenous case on the square lattice. In 
§4 we find the resultant generic formulas for the bond- 
percolation thresholds and compare our results to the 
very recent studies. In the last section, we conclude our 
study. 



II. CONVENTIONAL ANALYSIS 

Our analysis is based on the duality @, 0] , which is the 
simplest way to estimate the bond-percolation thresh- 
olds. We consider the bond-percolation thresholds for 
the lattice consisting of repetition of the unit cell as in 
Fig. [TJ Let us define p, r, s, and t as the inhomoge- 
nous probabilities to connect both ends of the assigned 
bonds. The conventional duality analysis can lead to 
the bond-percolation thresholds for the homogenous case 
p = r = s = t on the square lattice. In addition to the 
duality, the star-triangle transformation gives the bond- 
percolation thresholds for the inhomogenous case on the 
triangular and hexagonal lattices. First, let us review the 
conventional duality for convenience. 
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vice versa, after the dual transformation through the q- 
componcnt discrete Fourier transformation. When the 
dual lattice is the same as the original one, the system 
holds self-duality. For instance, the square lattice is the 
case. Then we can regard Z*(xq, x*, • • • ) as Z(xq, x*, • • • ) 
and can obtain the exact value of the critical point by the 
duality. We restrict ourselves to the case on the square 
lattice. Notice that -/Vb/2 = N$ on the square lattice. 
Let us extract the principal Boltzmann factors with edge 
spins parallel xq and Xq from both sides of Eq. © as 



FIG. 1: The unit cell of the inhomogenous bond-percolation 
problem on the square lattice. The unit cell is a part of the 
square lattice, which is covered with the dashed lines. The 
assigned values p, r, s and t are the probabilities to connect 
both ends on each bond on the unit cell. 

A. Duality 

We consider the g-state Potts model with the following 
Hamiltonian, 

H=-Y,Jii&(<t>i-<l>i), (!) 

m 

where Jij is the strength of interactions and takes differ- 
ent values as J p , J r , J s , and Jt, which will correspond to 
the probability assigned on the bonds. The summation 
is taken over all bonds, <5(x) is Kronecker's delta, and 4>i 
stands for the spin direction taking 0, 1, • • • , and q— 1. Let 
us estimate the critical point of the g-state Potts model, 
since it corresponds to the bond-percolation threshold in 
the limit of g -> 1 [IS]. 

We here assume the homogeneous case J = J p = 
J r = J s — J t . The duality exploits an inherent symme- 
try embedded in the partition function with the inverse 
temperature /3 as Z = H(ij) exp(/3 J<5(^ - (f>j)) — 
E* n<y)(l + vSfa - 0,)), where v = exp(/3J) - 1 0. 
Two different approaches to evaluate the partition func- 
tion, the low- and high-temperature expansions, can be 
related to each other by the g-component discrete Fourier 
transformation for the local part of the Boltzmann fac- 
tor, namely edge Boltzmann factor Xfc = I + v5(k) 
Specifically, each term in the low-temperature expan- 
sion can be expressed by Xk , while the high-temperature 
one is written in the dual edge Boltzmann factor x* = 
^ k Xfc exp(i2ivkl/q) / y/q. As a result, we obtain a double 
expression of the partition function by use of two differ- 
ent edge Boltzmann factors as 

Z(x ,x 1 ,---) = q Ns -^ L - 1 Z*(x* ,x* 1 ,---). (2) 

where Z* is the partition function on a dual lattice. Here 
Ns and Nb denote the numbers of sites and plaquettes, 
respectively. The unity in the power of q can be ignored 
in the following analysis. We obtain another system on 
the dual graph, on which each site on the original lat- 
tice exchanges with each plaquette on the dual one and 



(x ) Nb z( Ui ,u 2 , ■■■) = {x* f B z{ul,u* 2 , ■■■), (3) 

where z is the normalized partition function 
z(ui,U2,-") = Z/(x ) Nb and z(u\, u|, • • • ) = 
Zj (xq) Nb . We here define the relative Boltzmann factors 
u k = x fc /xo = 1/(1 + v) and u* k = x£/xq = v/(q + v). 
The well-known duality relation can be obtained by 
rewriting u* k in the same form as u k by use of v* as 
v/(q + v) = 1/(1 + v*), namely v* — q/v. Notice that 
the quantity v* has a different parameter K* from the 
original coupling K = f3J, which implies transformation 
of the temperature. We obtain the exact value of the 
critical temperature from the fixed point condition 
vf. = q under the assumption that a unique transition 
undergoes in the system. The limit q — > 1 can then 
give the bond-percolation threshold in the homogenous 
case p c = p = r = s = t on the square lattice through 
p c = v c /(l + v c ), namely p c = 1/2 [4, 5]. We can also 
derive the critical point by the following simple equality 

x =xS- (4) 

Indeed this equality gives v c = 1, namely p c = 1/2. 

For the case without self-duality, we can find an im- 
portant relation from v* = q/v. We can relate the prob- 
ability assigned on the bond on the original lattice to 
that on the dual one as p* = 1 — p in the limit q — > 1 
HQ. In other words, the probability p that both ends 
are connected on the original lattice is transformed into 
the disconnected probability on the dual lattice as I — p* 
and vise versa. We can rewrite this fact in terms of the 
relationship of the connectivity as 

P{AB) = P(A\B), (5) 

where the quantity on the left-hand side expresses the 
probability that A and B are connected, and that on the 
right-hand side stands for the probability that A and B 
are disconnected. The end points A and B in Fig. [^de- 
note the sites on the original lattice. On the other hand, 
A and B represent the sites on the dual lattice. Then 
the bond percolation threshold for the homogenous case 
on the square lattice can be represented by the following 
equality 

P(AB) = P(A\B). (6) 
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FIG. 2: Duality relation of the bond-percolation problem. 
The dotted line denotes the disconnected bond. The bold 
line represents the connected bond. The white circles denote 
the original sites. The black circles represent the dual sites 
(original plaquettes). 
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[Xq,xI,---) = Z{Xq ',x^ ',■■■) in Eq. ©. 
here use the renormalized-edge Boltzmann factor x k 
defined as 



We 

= (tr) 



*(tr) 



(7) 



where the product runs over i = p,r, and s for the three 
bonds on the unit cell of the hexagonal lattice, namely the 
down-pointing (up-pointing) star. We here assume the 
inhomogeneous system with vi — exp(/3Ji) — 1. We take 
the summation over the internal spin cf>o denoted by the 
black circle on the unit cell as in Fig. [3] The coefficient 
comes from that in front of the partition function 
on the right-hand side of Eq. ©. Notice that N$ is 
the same as the number of down-pointing (up-pointing) 
triangles N tr on the triangular lattice, and Nb = 3Ns- 
The subscript k denotes the configuration of the edge 
spins {4>i=p,r,s\ on the unit cell. On the other hand, 
we rewrite the original partition function in terms of the 
product of the edge Boltzmann factors as 



(tr) 
X k = 



- 4>o)vi) 



(8) 





FIG. 3: The inhomogeneous bond-percolation problem on the 
triangular and hexagonal lattices. The assigned values p, r 
and s are the connected probabilities assigned on each bond 
on the unit cell. On the hexagonal lattice, we put the dual 
probabilities p*, r* and s*, which are obtained after the dual 
transformation. The black circle on the hexagonal lattice rep- 
resents the internal site we sum over in the star-triangle trans- 
formation. 



B. Star-triangle transformation 

Let us consider the case on the triangular lattice. We 
here remove the homogeneous restriction that we impose 
above. We deal with the bond-percolation problem with 
the inhomogenous probabilities on the triangular lattice 
as depicted in Fig. [3] 

The dual transformation changes the triangular lat- 
tice into the hexagonal lattice. Then we cannot perform 
the same analysis as that in the case on the square lat- 
tice. We employ another technique to relate the hexago- 
nal lattice to the original triangular lattice. This can be 
achieved by the partial summation over internal spins at 
the down-pointing (up-pointing) star on the hexagonal 
lattice, namely star-triangle transformation Q . Then we 
can transform the partition function on the hexagonal 
lattice into that on another triangular lattice, namely 



The double expression of the partition function can be 
written as 



Z(xq \x[ ,■■■) — Z*(x Q 



y* i *(tr) *(tr) 



•)• (9) 



Similarly, let us extract the renormalized-principal Boltz- 



mann factors with edge spins parallel x^ 
both sides of Eq. © as 
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Notice that the number of the down-pointing (up- 
pointing) stars is the same as N tr . and is the 
normalized partition function z^ tr '(u^ tr ', u 2 tT \ ■ ■ ■ ) — 

Z/(x^) Nt ' and 2 ( tr )( U ; (tr) , M * (tr) ,...) = Z/(xf T) ) Ntt - 
We here define the renormalized-relative Boltzmann fac- 



tors u 



(tr) 



Xi, I Xq 



" and u^ tr ^ 



"(tr) 



/ * 



(tr) 



Simi- 



larly to the case on the square lattice, we put the simple 
equality as 



(.(tr) 



(11) 



This equality yields the critical manifold of the g-state 
Potts model as detailed in Appendix [5] By taking the 
limit q — > 1, we obtain the equality for the bond- 
percolation thresholds on the triangular lattice as 



where 



T(p,r,a) = 0, 



T(p, r, s) — prs — p — r — s + 1. 



(12) 



(13) 
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FIG. 4: The unit cell with four bonds for the inhomogenous 
case on the square lattice. The black circle denotes the inter- 
nal spins that we sum over, while the white ones are fixed as 
4>i = 0. 

If we perform the dual transformation on this equality, 
we find the solution for the bond-percolation thresholds 
on the hexagonal lattice as 

H{p\r\s*) = 1 (14) 

where 

H (p, r, s) = prs — rp — rs — ps + 1. (15) 

As shown above we can obtain the exact solution of the 
bond-percolation thresholds for several cases through the 
duality and the technique in conjunction with the star- 
triangle transformation. 

III. DUALITY WITH REAL SPACE 
RENORMALIZATION 

The duality with the star-triangle transformation, 
which is the partial summation of the unit cell on the 
hexagonal lattice, leads to the exact solution for the 
bond-percolation thresholds on the triangular and hexag- 
onal lattices as in (|T2"j) and (|14l) . Let us develop the sim- 
ilar analysis to the successful case on the triangular and 
hexagonal lattices. 

We start from Eq. ^ for the case on the square lattice. 
Notice that the edge Boltzmann factor is not enough to 
express the local property of the inhomogeneous system. 
Thus we consider to use the renormalized-edge Boltz- 
mann factor inspired by the star-triangle transformation. 
We take the square unit cell consisting of four bonds from 
both of the original and dual square lattices as in Fig. @] 
Let us take the product of the edge Boltzmann factors 
and perform the summation over the internal spin. The 



resultant quantity is written as 

4 sq) =E n (i+«<^-0o)). (i6) 

Similarly, we obtain the dual renormalized-edge Boltz- 
mann factor as 

We can then rewrite the relation obtained by the conven- 
tional duality ([2]) as 

Z(4 sq) ,x^,...) = ^ (sq) ,xt (sq) ,...). (18) 

We extract the renormalized-principal Boltzmann factors 
XQ Sq ' ) and with all edge spins on the unit cell parallel 

as 

(4 sq ^/v s ^(4 sq \4 sq) ,---) 

= ( x ;M ) ^/4,W( !i jW )U »(* 1 ) ) . . . ); (19) 

where the normalized partition function 

but z^)(4 sq) ,4 sq) ,-.-) = Z/(x^) Nb/4 and 
z« (sq) ,^ (sq) ,---) - Z/(x* (sq) ) N */\ We here de- 
fine the renormalized-relative Boltzmann factors 
u («0 = 4 sq »/4 sq ) an d u ^ sq ) = xl (sq) /x* (sq) . Then we 
impose the following equation to identify the location of 
the critical point 



The direct evaluation of this equality in the leading order 
of e where q = 1 + e gives the formula for the bond- 
percolation thresholds, as detailed in Appendix iBl 

Y[(l + Vi)C(p,r,s,t) = 0. (21) 

i 

where 

C(p,r, s,t) = 1 — pr — ps — rs — pt — rt — st 
+prs + prt + rst + pst. (22) 

It is reasonable that a unique transition undergoes if we 
tune the temperature for the inhomogeneous interactions 
as J p , J r , J s , and Jt, which correspond to the probabili- 
ties of the bond-percolation problem in the limit q 1. 
Therefore the singularity of the free energy should be 
unique for change of the temperature. The duality can 
then identify the location of the critical point. Therefore 
we conclude that C(p,r,s,t) — gives the exact bond- 
percolation thresholds for the inhomogeneous case on the 
square lattice. 

The equality C(p, r, s, t) = was originally conjectured 
@, confirmed numerically with high precision and de- 
rived in a different way [9]. The proof of validity of 
Eq. (j2"2"j) has been very recently established 0. It is 
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not simple to show the validity of Eq. (|22|) . since it 
is based on the indirect analysis via considerations of 
the bond-percolation problem on different lattices. The 
present analysis demonstrates the more straightforward 
analysis. Without recourse to the duality, the real-space 
rcnormalization group analysis can give the exact bond- 
percolation thresholds for the homogenous case but fails 
into the approximations for the inhomogenous case [1? 
]. By virtue of the duality, we can here find the exact 
answer for the critical point while we stand on the fixed 
point of the renormalization group. 

The duality with real-space renormalization as shown 
above is essentially the same as the profound technique 
in the analysis of the random spin system, in particular 
spin glasses [ll|, ■ For several models in the random 
spin system, the dual transformation cannot relate the 
original system to the same one with a different tem- 
perature as the case for the g-state Potts model despite 
existence of self-duality of the lattice. In these cases, we 
recover the self-duality of the random spin models via 
real-space renormalization over larger range beyond the 
unit cell, namely summation over several internal spins 
by taking larger size of the cluster in order to find the 
correct fixed point in a relatively wide space of param- 
eters as well as the temperature. Then we impose the 
following condition to estimate the critical point simi- 
larly to the conventional analysis by the duality ((4]), its 
combination with the star-triangle transformation (ITT1) . 
and the above analysis ((20] 



(23) 
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where iEq and 



are renormalized principal and dual 
Boltzmann factors on the cluster. The size of the cluster 
is denoted by b (b = means the simple duality with- 
out renormalization). It is examined that, if b is taken 
to be a large value, we can recover the self-duality fol- 
lowing the concept of renormalization [ID, [l4| , and Eq. 
(I23p can give a precise estimation of the critical point 
11]. When we deal with the random spin system on 
the square lattice, the estimation by setting b — 1 (four- 
bond cluster as depicted in Fig. [5]) often attains satis- 
fiable precision. If one wishes to enhance the precision, 
the systematic improvement is achievable by increase of 
b (i. e. 6 = 2 16-bond cluster as in Fig, Indeed we 
can also estimate the bond-percolation thresholds from 
the critical points via analysis on the bond-dilution Ising 
model, which is a typical model in the random spin sys- 
tem. In Appendix \C\ we demonstrate the rederivation of 
the formula C(p,r, s,t) = through the duality analy- 
sis with real-space renormalization for the bond-dilution 
Ising model with inhomogenous distribution. Equation 
(f2"3")) yields the formula C(p,r,s,t) = for the bond- 
percolation threshold for the inhomogeneous case on the 
square lattice. Below we examine the obtained result 
through the duality with real-space renormalization in 
several points of view. Beyond the case on the square 
lattice, we try to apply several similar systems. 
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FIG. 5: The cluster for the duality with real-space renormal- 
ization on the square lattice. The white circles denote the 
fixed spins for evaluation of the principal Boltzmann factors. 
The black circles express the spins we trace over (similarly 
to the star-triangle transformation) to obtain the principal 
Bolztmann factors. 



IV. GENERIC FORMULAS ON BOND 
PERCOLATION THRESHOLDS 

When we analyze the homogenous case on the square 
lattice, the formula of the bond-percolation threshold is 
a relation on the unit cell with the two-terminal struc- 
ture (a single bond) as in Eq. (|6]). On the other hand, 
the unit cell where we perform the analysis consists of 
the three-terminal structure for the case on the triangu- 
lar and hexagonal lattice as depicted in Fig. [3] As in 
this case, for the lattice consisting of repetition of the 
three-terminal unit cell, the generic formula of the bond- 
percolation thresholds is known to be [15j ■ 



P(ABC) = P(A\B\C), 



(24) 



where the quantity on the left-hand side expresses the 
probability that the end points A, B, and C on the three- 
terminal unit cell are all connected, and that on the right- 
hand side stands for the probability that none of A, B, 
and C are connected. 

We obtain the critical manifold (fT2")) for the bond- 
percolation thresholds on the triangular lattice from the 
above formula ([24]) . We here demonstrate the reduction 
to Eq. (p~2]) from Eq. (|24|) . In the case on the triangular 
lattice, let us write down all terms included in P(ABC) 
as 

P(ABC) = prs+pr(l - s) + p(l - r)s + (1 -p)rs 
= pr + sp + rs — 2prs. (25) 

On the other hand, the probability that none of the end 
points are connected is 

P(A\B\C) = (i- p )(i_r)(l-a) 

= 1 — p — r — s + pr + rs + sp — prs. 

(26) 

Thus we can obtain Eq. (JTSJ) from Eq. (|24|) . 

In addition, on the hexagonal lattice, Eq. (|24|) can 
be reduced to Eq. (|14l) . On the hexagonal lattice, the 
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FIG. 6: Transformations into the hexagonal and triangular 
lattices from the square lattice. 

left-hand side of Eq. ([24| is written as 

P(ABC) = prs. 

The right-hand side can be given as 

P{A\B\C) = (i-p)(i_ r )(l-«)+p(l-r)(l-a) 
+ (1 - p)r(l - s) + (1 -p)(l - r)s 
= 1 — pr — rs — sp + 2prs. (27) 

Equation (J23]) reproduces Eq. ([T3]). 

We here give a generic formula for several lattices con- 
sisting of repetition of the four-terminal unit cell as in 
Fig. |5] (left). The analysis as detailed in Appendix [C] 
provides the generic formula for the four-terminal unit 
cell as 

P(ABCD) + P(BCD\A) + P(ACD\B) 

+P{ABD\C) + P{ABC\D) = P(A\B\C\D), 

(28) 

where P(BCD\A) is the probability that BCD connects 
with each other while A is disconnected, and the other 
quantities follow the same manner. We can reproduce 
Eq. (p4|) by reduction to the three-terminal unit cell (in 
particular hexagonal lattice) by removing the single bond 
from four bonds as in Fig. [5] (left). It means that all 
connected probabilities to D vanish as P(ABCD) = 
since D can not be connected, and we omit depen- 
dence on D for the disconnected probability with D as 
P{ABC\D) = P(ABC). 

For the case of the four-terminal unit cell, we can 
obtain another generic formula. By the conventional 
duality, we relate the bond-percolation problem on the 
original square lattice to that on the dual square lat- 
tice through the duality relation. The probabilities ex- 
pressing connectivity of the edge sites are then changed 
as P(ABCD) = ' P{A\B\C\D) and P(D\ABC) = 
P{DC\A\B) similarly to Eq. ©. Another generic for- 
mula for the four-terminal unit cell as in Fig. (|6]) (right) 
can be expressed as 

P(A\B\C\D) + P(AB\C\D) + P(BC\D\A) 
+P(CD\A\B) + P(DA\B\C) = P(ABCD), 

(29) 

which is detailed in Appendix [Cj Here let us again reduce 
the above equality to the case of the three-terminal unit 




FIG. 7: Covering by the two four-terminal unit cells. The 
shaded squares express the unit cells. The left panel denotes 
the original square lattice, and the right one depicts the dual 
lattice. 

cell. This can be achieved by eliminating the terms asso- 
ciated with the disconnected probabilities with D since 
they are always connected as in Fig. |6] (right). In ad- 
dition we omit the dependence on D. Then Eq. (12"9")) 
recovers P(ABC) = P(A\B\C), namely Eq. dHJ). The 
difference between the four-terminal unit cells associated 
with Eqs. (f2"8"| and (|2"5|) comes from the tiling manner to 
cover the whole lattice. The former case is the full tiling 
of the unit cell. On the other hand, the latter case is the 
checker-board tiling as in Fig. [7J 

In order to show the efficiency of the above generic 
formula, we take a fascinating instance of the applica- 
tion beyond the case of the square lattice. By the above 
generic formula (f2"9"| . we can recover the equality for the 
bond-percolation thresholds on the bow-tie lattice, which 
is dealt with to prove the validity of Eq. (J24J) 7]. The 
four-terminal unit cell of the bow-tie lattice is shown in 
Fig. [5J Equation (|2"9"]l can be then reduced to 

u{prs(\ — t) + pr(\ — s)t + p(l — r)st + (1 — p)rst 
+P(l - r)(l - s)t+p(l - r)s(l - t) 
+ (1 - p)r(l - s)t + (1 - p)rs(l - t) 
+prst} + (1 - u)C(p, r, s, t) 
= C(p, r, s, t) — u(l — pr — st — prst) = 0. (30) 

This equality has been given by combination of the re- 
sults for the three-terminal unit cells by splitting of 
the four-terminal unit cell to two triangles as in Refs. 
0, B Ell ■ Then the combination of the duality and star- 
triangle transformation yields the exact solution of the 
bond-percolation thresholds on the bow-tie lattice. 

As another interesting but wrong instance of applica- 
tions, let us apply the generic formula (f2T)| to the bond- 
percolation problem on the Kagome lattice by consid- 
ering the four-terminal unit cell with up-pointing and 
down-pointing triangles as in Fig. [9l Here we take the 
homogenous case for simplicity. We then obtain the fol- 
lowing polynomial from the generic formula (|29[) as 

1 - 3p 2 - 6p 3 + 12/ - 6p 5 + p 6 = 0. (31) 

The solution is p c — 0.52442971, which is known to be 
an approximate estimation (Wu's conjecture Q), while 
a numerical evaluation gives p c = 0.52440502(5) p7| . 



7 



_ 




FIG. 8: The bow-tie lattice. The shaded squares express the 
four-terminal unit cells. 




D C 



FIG. 9: The Kagome lattice. The four-terminal unit cell con- 
sists of the up-pointing and down-pointing triangles. 



The reason why the above equality yields an wrong value 
comes from the solvability by the duality and the stat- 
triangle transformation. Our generic formula (1291) then 
fails to give the exact answer but possibly an approxi- 
mate estimation when the system lacks the solvability. 

However our formulation by the duality with real- 
space renormalization can give a more precise value of 
the bond-percolation threshold on the Kagome lattice. 
The above case is the same as that we used to find in 
study on random spin systems [U EH • The generic for- 
mulas as in (|28p and (|29p can be given by Eq. (|23p for 
6 = 1 through analyses of the random spin systems as de- 
tailed in Appendix [Cj We can obtain the precise value of 
the bond-percolation threshold by considering the larger 
cluster over the four-terminal unit cell, namely 16-bond 
cluster (b = 2) and more. Although we here cease the dis- 
cussion on the bond-percolation problem on the Kagome 
lattice since its threshold is out of scope in the present 
study, we remark that several recent development on 
this issue. The similar analysis is proposed in context 
of the graph polynomial as demonstrated in Refs. (]~8l — 
[20j . The idea is essentially based on consideration on the 
large cluster over the four-terminal unit cell to estimate 
the bond-percolation thresholds on the Kagome lattice. 
The method has indeed succeeded in giving very pre- 
cise estimations of the bond-percolation threshold on the 
Kagome lattice as p c = 0.52440500(1) These fact 

suggest that two independent methods developed in ran- 
dom spin systems and graph polynomial would be closely 
related to each other through the bond-percolation prob- 



lem. We hope that the future study reveals more clear 
relationship between their different realms. 



V. CONCLUSION 

In the present study, we rederived the exact solution 
for the bond-percolation thresholds in the inhomogenous 
case on the square lattice by use of the duality with real- 
space renormalization, which is a generalized analysis of 
the star-triangle transformation. In addition, we obtain 
two different generic formulas depending on the tiling 
manner of the four-terminal unit cell to cover the whole 
lattice. Both equalities can be reduced to the known 
formula for the three-terminal unit cell, which includes 
the triangular and hexagonal lattices. The application 
of the generic formula reproduces the exact solution on 
the bow-tie lattice and the well-known approximate solu- 
tion on the Kagome lattice. The further analysis possibly 
gives more precise value of the bond-percolation thresh- 
old on the Kagome lattice. 

The duality analysis with real-space renormalization 
shown in the present study is essentially the same as the 
special technique, which has developed in context of spin 
glass theory. The method has been useful to describe 
the precise phase boundary [ll|, EH IU- The straight- 
forward rederivation of the existing results on the bond- 
percolation problem are given in the different context im- 
plies existence of a fascinating theoretical connection be- 
tween different realms, graph polynomial and theory of 
spin glasses. 

We emphasize high nontriviality of our results shown 
in the present study. The exact solutions for finite di- 
mensional many-body systems have been rare in spite 
of the long-year efforts. However the situation begins 
to change by development of the duality analysis, which 
is found to be applicable to a relatively broad class of 
problems, namely spin glasses and inhomogenous perco- 
lation problems. We hope that the duality analysis with 
real-space renormalization would play an essential roll to 
understand the nature of the many-body systems as the 
conventional duality proposed by Kramers and Wannier 
contributed to establishment of the Onsager solution [22j ■ 
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Appendix A: Derivation of Eq. (|12[) from Eq. 

We evaluate Eq. (fTTj) in this appendix. From the def- 
inition, we write down Eq. (1111) as 

ri( i + ^ = \ { - !) n * + U(i + } ■ ( ai ) 

This is the critical manifold of the g-state Potts model on 
the triangular lattice. Let us take the leading term of e 
of q = 1 + e for obtaining the bond-percolation thresholds 
on the triangular lattice. 



*(tr) _ *(tr) 



-«nc+«.)(-'+Ern;+ni^)- 

(A2) 

By rewriting each in terms of the probability assigned 
on each bond as p = + u p ) etc., we reach 



JJ(l + « 4 )T(p,r,8) = 0, 

i 

which is reduced to Eq. (fT2|). 



(A3) 



Appendix B: Derivation of Eq. (|22[) 

We here demonstrate the detailed evaluation of Eq. 
(El) from Eq. (I2TH) . We can write the difference between 
the left and right-hand sides of Eq. (j2"U)) by use of defi- 
nition of the renormalized-edge Boltzmann factors as in 
Eqs. (fill) and (ED as 

*( s q) „(sq) 



^{•- 1+ nH)} 
-|«-i+n(i+«i)|. 



(Bl) 



We take the leading term of e of q = 1 + e. In advance, 
we evaluate the following quantities 

*H«- 1+ nK)} 

-naWi-fc+.E^IIifs) 

i \ i i / 



(B2) 



and 



jg-i+nu+zoj 
=nc i +««)( i + e iiT^) ( B3 ) 



Therefore Eq. (IB1[) can be reduced to 

n( i+ ^) 

^{-^i:T^- + (n».- 1 )nT^-} 

= JJ(l + «i)C(p,r,a,i) 

i 

We reproduce Eq, 



(B4) 



Appendix C: Alternative way to Eq. (|22[) 

We show an alternative way to give Eq. (f2"2"j) with 
recourse to the bond-dilution Ising model on the square 
lattice. We consider the following Hamiltonian 



H— jjjSjSj, 

(ij) 



(CI) 



where Si stands for the Ising spin taking ±1, and 
stands for the random coupling following the distribution 
functions 



P P (J,) = p5(J iS - J) + {l-p)5{J l3 ). 



(C2) 



We also define P s , Pt and P u for each bond. 

In random spin system, we need to take the configu- 
rational average of Jij to evaluate the free energy. We 
often employ the replica method to perform the configu- 
rational average. Instead of the averaged logarithm of the 
partition function (free energy), we analyze the averaged 
power following the well-known identity as 

[Z n ] - 1 



[log Z] = lim 

n— >0 



(C3) 



where [• • • ] expresses the configurational average. Ini- 
tially we deal with the replicated system by setting n as 
a natural number. At the final step of analysis, we take 
the limit of n — > 0. We then regard the averaged power of 
the partition function [Z n ] as the effective partition func- 
tion written as Z n (the replicated partition function). 

Let us perform the duality analysis with real-space 
renormalization by dealing with the effective partition 
function. The effective partition function consists of the 
following edge Boltzmann factor as 



[J exp(K nj S?S?) 



a = l 



(C4) 



where K = (3 J, and Ty takes or 1 and expresses the 
existence of the interaction. The superscript a runs from 
1 to n standing for the index of the replicas. On the other 
hand the dual edge Boltzmann factor is defined as 



1 

.3 



(C5) 
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Let us take the cluster with four bonds as in Fig. [5] 
to evaluate Eq. (f2"3f for 6=1. In order to evaluate 
the renormalized-edge Boltzmann factors, we fix the edge 
spins to Si = 1 on the cluster and sum over the internal 
spin similarly to the star-triangle transformation. The 
renormalized-principal Boltzmann factor is written as 



T (l) _ 

x — 



^[]cxp(^r 4 5 ) 



So i 



(C6) 



where the product runs over i = p,r,s and t. The dual 
renormalized-principal Boltzmann factor is given by 



*(i) 



S i 



(C7) 



In order to identify the location of the bond-percolation 
thresholds, we consider K — > oo. We obtain 



log (V^-rha-.) (i + l[r^j 



2 log 2. 



(C9) 



First, let us take the homogenous case p = r = s = t 
Equation (|CJ9[) becomes 



P A - V(l -p) + 4p(l - pf + (1 - p) = 0. (CIO) 



This implies 



Taking n — > in Eq. (|23|) . we obtain the following for- 
mula 



= 2 log 2. 

(C8) 



p 4 + V(l -p)= 4p(l - pf + (1 - p) 4 . (Cll) 



Let us obtain the generic formula for the inhomogenous 
case (HU). From Eq. (EH) we find 



prst + {prs(l — t) + pr(l — s)t + p(l — r)st + (1 — p)rst} 

- {p(l - r)(l - «)(1 - t) + (1 -p)r(l - - t) + (1 - r)a(l - t) + (1 -p)(l - r)(l - a)t} 

-(l-p)(l-r)(l-s)(l-t) = 0. (C12) 



By simplifying the above equality, we reproduce Eq. (|22[). 

We can find the general formula (|28p as follows. The first term in Eq. (|C12|) prst corresponds to P(ABCD) in Eq. 
([28]) . The following four terms in the first line is P(ABC\D), P{ACD\B), P (ABD \C), and P(BCD\A), respectively. 
The remaining terms become —P(A\B\C\D). Therefore the above equality (|C12[) suggests the generic formula (|28|) . 

On the other hand, the simple duality as p* = 1 — p reduces Eq. (IC12[) to 

p*r*s*t* + {p*r*s*(l - t*) + p*r*{l - s*)t* +p*(l - r*)s*t* + (1 - p*)r*s*t*} 

- { P *(l - r*)(l - s*)(l - t*) + (1 - p*)r*(l - s*)(l - t*) + (1 - - r*)s*(l - t*) + (1 - p*)(l - r*)(l - s*)f} 
-(1 - p*)(l - r*)(l - s*)(l - t*) = 0. (C13) 

I 

Then the collection of all the terms in the first line be- and — P(CD\A\B), respectively. The last term is noth- 
comes P(ABCD). Each term in the second line corre- ing but — P(A\B\C\D). We thus obtain another generic 
sponds to -P(AB\C\D), -P(BC\D\A), -P(DA\B\C), formula flgj). 
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